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It is shown that the generalized discrete nonlinear Schrodinger equation can be reduced in a 
small amplitude approximation to the KdV, mKdV, KdV(2) or the fifth-order KdV equations, 
depending on values of the parameters. In dispersionless limit these equations lead to wave breaking 
phenomenon for general enough initial conditions, and, after taking into account small dispersion 
effects, result in formation of dissipationless shock waves. The Whitham theory of modulations of 
nonlinear waves is used for analytical description of such waves. 

PACS numbers: 05.45.Yv, 05.90.+m 

I. INTRODUCTION 

Dissipationless shock waves have been experimentally observed or their existence has been theoretically predicted 
in various nonlinear media - water pj , plasma 0] , optical fibers Q , lattices . In contrast to usual dissipative shocks 
where combined action of nonlinear and dissipation effects leads to sharp jumps of the wave intensity, accompanied by 
abrupt changes of other wave characteristics, in dissipationless shocks the viscosity effect is negligibly small compared 
with the dispersive one, and, instead of intensity jumps, the combined action of nonlinear and dispersion effects leads 
to formation oscillatory wave region (for review see e.g. Since intrinsic discreteness of a solid state system gives 
origin to strong dispersion, which can dominate dissipative effects in wave phenomena, it is of considerable interest 
to investigate details of formation of dissipationless shocks in lattices. 

As a model, in the present paper we choose the general discrete nonlinear Schrodinger (GDNLS) equation 

iq„ + (1 - r)\q„\ 2 )(q n +l + q n -l - 2q n ) + 2(p 2 - \q n \ 2 )q n = (1) 

introduced by Salerno 0, Q ■ Eq. Q turned out to be an important model not only because of its property to provide 
a one-parametric transition between an integrable Ablowitz-Ladik (AL) model and the so-called discrete nonlinear 
Schrodinger (DNLS) equation (r) = 1 and r\ — 0, respectively), but also because of a number of physical application 
for a review of which we refer to |8(- 

Obviously, Q has a constant amplitude solution q n = p. It was shown in [flLllCj. that in a small amplitude, \a\ <C p, 
and long wave (so that the discrete site index n can be replaced by a continuous coordinate x) limit, evolution of 
small amplitude perturbations with respect to this constant background, 

q n (t) = (p + a(x, t)) exp(-i0(a;, t)), (2) 

is governed by the Korteweg-de Vries (KdV) equation for the amplitude a(x, t): 



2(3 - y 7 !^ ? 
— -aa x H 



at , aa x H — [3(1 - r\p ) - p \a xxx = (3) 



which is written in the reference system moving with velocity 2p\J\ — r/p 2 of linear waves in dispersionless limit. It 
is well known (see, e.g. [l|) that if the initial pulse is strong enough, so that the nonlinear term dominates over the 
dispersive one at the initial stage of the pulse evolution, then the dissipationless shock wave develops after the wave 
breaking point. The theory of such waves, described by the KdV equation, is well developed (see, e.g. @)- Existence 
of the respective shock waves for model has been predicted analytically and observed in numerical simulations in 
0,^3- Moreover, as it is shown in ^lj a nonlinear Schrodinger equation of a rather general type can bear "KdV- type" 
shock waves. In this context the results presented below in the present paper although being mostly related to model 
Q display some general characteristic features of a lattice of a nonlinear Schrodinger type. 

The coefficients of (jSJ depend on two parameters rj and p and can vanish at special choice of these parameters, 
so the KdV equation looses its applicability for these values of r\ and p. Physically acceptable values of r\ and p are 
limited by the inequalities 

< T) < min{l, l/p 2 }, 0<p<oo, (4) 
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FIG. 1: Diagram showing different continuous limits of the GDNLS equation. The abbreviations are explained in the text. 



and this region of the (p, 77) plane is depicted in Fig. 1. 
Along the line 



the coefficient before the nonlinear term in @ vanishes. This means that implied in derivation of the KdV equation 
approximation limited to only quadratic nonlinearitics fails and higher nonlinearities have to be taken into account 
for accurate description of the wave dynamics. Hence, one may expect that the modified KdV (mKdV) equation with 
cubic nonlinearity arises for values of p and r\ related by (J5J • 
Along the line 
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V = ^ ~ o (6) 



the coefficient before the dispersion term in vanishes what means that higher order dispersion effects have to be 
taken into account. In this case one may expect that evolution equation for a(x,t) contains quadratic nonlinear term 
and linear dispersion term with fifth-order space derivative of a(x,t). 
The most interesting point corresponds to the values 

r) = 1, p = V3/2 (7) 

when both nonlinear and dispersion coefficients vanish. Since Eq. Q with rj — 1 coincides with the completely 
integrable Ablowitz-Ladik equation, one may expect that in a small amplitude approximation this equation with the 
parameters equal to reduces again to a completely integrable equation. The KdV equation Q with 77 = 1 is valid 
for all interval < p < 1 except for some vicinity of the point p = Vo/2, and therefore one can suppose that at the 
point one has to obtain the second equation of the KdV hierarchy KdV(2) in which the higher order nonlinear 
and dispersion effects play the dominant role. 

The aim of this paper is two-fold. At first, in Sec. [H] we derive the evolution equations for the whole region 0} and 
show that along the line © the small-amplitude approximation reduces to the mKdV equation fSec III B|) . along the 
line © to a nonlinear equation with the dispersion of the fifth-order fSec III C]) . and at the point to the KdV(2) 
equation fSec III Dft . The last result sheds some new light on the nature of higher equations of the KdV hierarchy — 
they arise as small amplitude approximations to completely integrable equations, if lower orders of nonlinear and 
dispersion contributions vanish at some values of the parameters of the equation under consideration. 

The second aim of the paper is to develop a theory of shock waves fSec llllf) for the mKdV fSec lIII All and KdV(2) 
fSec lIIIBjl equations analogous to that developed earlier for the KdV equation. This theory permits one to described 
in details the behavior of shocks after the wave breaking point for different values of the parameters entering into 
Eq. ft). The results are summarized in Conclusion. 
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II. SMALL AMPLITUDE APPROXIMATION 

Using ansatz and replacing a discrete index n by a continuous variable x, we rewrite Eq. Q as 

ia t + (p + a)(j)t — 4:p 2 a — 6pa 2 + [1 - ri(p + a) 2 ] 
x{(p + a(x + 1, t)) exp[-i(>(a: + 1, t) - (f>(x, t))} 
+ (p + a(x - 1, t)) exp[-i((/>(x - l,t) - 4>{x, t))]} = 0. 

In the linear approximation, this equation yields for the harmonic wave solution 

a(x,t) cx ex.p[i(Kx — fit)], ^(x,t) cx exp[i(Xa; — fit)], 

the dispersion relation 0, 



(8) 



fi = ±4-\/T - ??P 2 sin 



p 2 + (1 - ?7p 2 )sin 2 



if 



1/2 



= ±2pVl -r]p 2 K- 



3(1 - r^ 2 ) - p 2 
24p 2 



2 + 0(K 4 ) 



(9) 



where expansion in powers of K corresponds to taking into account different orders of the dispersion effects. In the 
lowest order, when the dispersion effects are neglected, linear waves propagate with constant velocity 



±2py/\ - Vp 2 . 



(10) 



To evaluate contribution of small (for \a\ <C p) nonlinear effects, it is convenient to introduce a small parameter e ~ a 
and pass to such scaled variables in which nonlinear and dispersion effects make contributions of the same order of 
magnitude into evolution of the wave. Since the choice of these scaled variables depends on values of the parameters, 
we shall consider the relevant cases separately. 



A. KdV equation 

For the sake of completeness we start by reproducing briefly some results of 0. We expand a(x± 1, t) and <fi(x± 1, t) 
into the Taylor series around x, introduce scaling indexes 

a~£, t~£" a , x~ £ - , </>~e 7 , (11) 

and demand that in the reference frame moving with velocity <|10|) of linear waves the lowest quadratic nonlinearity 
has the same order of magnitude as the second term in the expansion of the dispersion relation, 

a ~ <Px: <H ~ o-a x ~ a xxx , 

which yield 

a 7j, (3 'f ^ • 

Thus, the scaled variables have the form 

r = e 3/2 t, £ = e 1/2 (x + vt), v = ±2py/l - rjp 2 , (12) 
and a(x, t) and (f)(x, t) should be looked for in the form of expansions 



a = ea^+e 2 a^+e 3 a^ + ..., 

= e 1 /20(l) +£ 3/2^(2) +e 5/2^(2)_ _ 

Then in the lowest order in expansion of Eq. JSJ in powers of e we obtain the relationship 



6« = 
5 v 



(13) 



a (1) , (14) 
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where we have chosen the upper sign of v in l|12fl , and in the next order the KdV equation © written in terms of the 
scaled variables. The nonlinear term changes its sign at (J3J and the dispersion term changes its sign at 10- Hence, 
there can be as bright solitons against a background (a(x, t) > 0) or dark solitons (a(x, t) < 0) of the GDNLS equation 
approximated by the KdV equation ©. In both cases it is possible to pass to new dependent variable 

--T^^[3(l- W 2 )-P> (15) 



and change time as t — > ^\ 2 ^, P [3(1 — f]p 2 ) — p 2 ] t such that the KdV equation takes the standard form 

u t + 6uu x + u xxx = 0. (16) 

B. mKdV equation 

At r\ = 3/4p 2 the quadratic nonlinearity in Eq. @ becomes zero and in order to describe the wave in the vicinity 
of the breaking point the scaling should be chosen to take into account cubic nonlinearity which is now must have the 
same order of magnitude as a xxx , 

a ~ (f> x , at ~ a 2 a x ~ a xxx . 

Then, using scaling we find 

a = 3, = 1, 7 = 0, 
so that instead of l|12(l we have the following scaled variables 

T = e 3 t, £ = e(x + vt), v = ±p, (17) 
where value of velocity is found by substitution of JSJ into (|1U|) , and instead of the second expansion in (|13f) we have 

= 0(l) +e(/) (2) +£ 2^(3) + __ 

Then in the lowest order in expansion of Eq. JSJ) in powers of e we obtain 

# 5 = -a« (18) 

what coincides with Eq. 114|) after substitution of v — p. In the next order we get the relationship 

=4a ( 2) + 6 a(1)2; (1Q) 
4 P 

and finally in the highest relevant order we obtain the mKdV equation 



<" + { i + 7) ""^ + W P { 1 ~ iJ P) * = a (20) 

Since along the line © we have p > V3/2, the coefficient before the dispersion term is always negative and hence 
Eq. I|2()[| can be transformed to the following standard form 



u t - 6u 2 u x + u xxx = 0, (21) 



where 



and wc rcnormalizc the time variable 



gV' p^Lu (22) 
4V(V + 21) V ' 



(23) 
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C. Fifth-order KdV equation 

At 77 = 1/p 2 — 1/3 the first order dispersion effects in Eq. disappears and in this case scaling should be chosen 
so that the quadratic nonlinearity has the order of magnitude of a^ v \ 



which yields 



(V) 

at ~ aa x ~ a , 



a=|, 0=1 7-^ 



45 k 41 / 45 
so that the scaled variables are given by 

T = e 5 /\ t = e x '\x + vt), v = ±^=, (24) 

v3 

where velocity v is found by substitution of JSJ into (| XOp . and now the condition of cancellation of terms in the second 
order demands that expansions of a and (/> have the form 

a = sa^ +e 3 / 2 a^ +e 2 a^ + ..., 

= e 3 /y 1} + £ 5 /y 2) + e ? /v 2 > . . . (25) 

Then in the lowest order in expansion of Eq. JSJ in powers of e 1 / 2 we get again 

= V3 a(1) 
5 P 

what can be obtained from Eq. I|14|l by substitution of ©. In the next order we get the relationship 



( 2 ) 2\/3 (2) 1 (i 

Pi = a v ; r=~ a ff 



« (2) - -^4f (27) 



and finally in the highest relevant order we obtain the equation 

Here the nonlinear term can be obtained from the corresponding term in the KdV equation © by substitution of © 
and dispersion term reproduces the expansion of the dispersion relation at the same value of 77. 



D. KdV(2) equation 

At the point Q both the quadratic nonlinear and first order dispersion terms disappear, so that now the cubic 
nonlinearity must be of the same order of magnitude as aS v \ 

a~<^) x , a t ~ a 2 a x ~ a^ v \ 

which yield 



2 > @ T 2 ' 



and the scaled variables are given by 

T = e 5 / 2 t, Z = s 1 / 2 {x + vt), w = ±^, (29) 

where tj is the velocity of linear waves at the point Q. The variables a(x,t) and <p(x,t) have the same form of 
expansions l|13fl as in the KdV equation case. In the lowest order we obtain, as one should expect, Eq. (|18|) : in the 
next order we get the relationship 

0« =4a< 2 3+4V3a« a -§ag>; (30) 
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and at the last relevant order we obtain the equation 

a« + 16V3a^ - I - § + « = 0. (31) 

As one should expect, the main nonlinear term here coincides with that of the mKdV equation H20JI at the point J2J, 
and linear dispersion term with the corresponding term of the fifth-order KdV equation (|28[l at the same point. 
By means of replacements 

a (1) = \=u, r = -30V3i, £ = x 

8^3 

Eq. Ij31|l can be transformed to standard form of the second equation of the KdV hierarchy — KdV(2) (see, e.g. Q): 

u t = ^u 2 u x + bu x u xx + \uu xxx + \u (v) (32) 

Thus, the completely integrable AL equation reduces in the small amplitude approximation either to the KdV 
equation beyond some vicinity of the point J7J , or to the second equation of the KdV hierarchy at the point Q , so 
that approximate equations remain completely integrable in both cases. This observation suggests that the property 
of complete integrability preserves in framework of singular perturbation scheme, which is a known phenomenon (see 
e.g. 13]). Then higher equations of some hierarchy may arise as approximate equations. This happens if at some 
values of the parameters of the underline completely integrable problem nonlincarity and, hence, dispersion of lower 
equations of the hierarchy vanish. This phenomenon can be viewed as physical meaning of the higher equations of 
hierarchies of integrable equations. 



III. DISSIPATIONLESS SHOCK WAVES 



In dispersionless limit when dispersion effects can be neglected compared with nonlinear ones, all derived above 
equations reduce in the leading approximation to the Hopf-like equation 

u t + u n u x = 0, (33) 

where n = 1 for the KdV and fifth-order KdV equations and n — 2 for mKdV and KdV(2) equations. It is well-known 
(see, e.g. 0) that Eq. Q with general enough initial condition leads to formation of the wave breaking point after 
which the solution becomes multi- valued function of x. This means that near the wave breaking point one cannot 
neglect the dispersion effects. If we take them into account, then the multi- valued region is replaced by the oscillatory 
region of the solution of the full equation. This oscillatory region is called dissipationless shock wave and its analytical 
description is the aim of this section. 

Existing theory of dissipationless shock waves can be applied in principle to completely integrable equations only. 
Among equations derived in the preceding section, however, fifth-order KdV equation (12811 docs not belong to this 
class. Fortunately, just this case of zero first-order dispersion was studied numerically in [QL Ilflj . We also bear in mind 
that the dissipationless shock waves of the KdV equation are already described in literature (see e.g. @). Therefore 
we shall not consider this equation here and concentrate our attention on the completely integrable models. 

The analytical approach is based on the idea that the oscillatory region of the dissipationless shock wave can be 
represented as a modulated periodic solution of the equation under consideration. If the parameters defining the 
solution change little on a distance of one wavelength and during the time of order of one period, one can distinguish 
two scales of time in this problem - fast oscillations of the wave and slow change of the parameters of the wave. Then 
equations which govern a slow evolution of the parameters can be averaged over fast oscillations what leads to the 
so-called Whitham equations y| and their solution subject to appropriate initial and boundary conditions describes 
evolution of the dissipationless shock wave. This approach was suggested in ^1 an d now it is well developed for the 
KdV equation case (see, e.g. |j|). The results of this theory can be applied to shocks in GDNLS equation when it is 
reduced to the KdV equation ® or P^l* ■ Since this theory is presented in detail in , we shall develop first analogous 
theory for the mKdV and KdV(2) equations and after that compare the results obtained for different equations. 



A. Dissipationless shock wave in the mKdV equation l|21|l 

At first we have to express a periodic solution of the mKdV equation (|21|l in a form suitable for the Whitham 
modulation theory. Such a form is provided automatically by the finite-gap integration method which is used here to 
find the one-phase periodic solution of the mKdV equation. 
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The finite-gap integration method (see, e.g. Q ) is based on the complete integrability of the mKdV equation, 
following from a possibility to represent this equation as a compatibility condition ty xt = ^> tx of two linear systems 



^ x = U*, * t = V*, * = 



U 2 r 



|TJ= iAiu y= M (34) 



-in iX J ' 1 C 

A = — 4iA 3 — 2m 2 A, 5 = 4iuA 2 — 2u x X — iu xx + 2iu 3 , C = — 4iuA 2 — 2u x X + iu xx — 2iu 3 

where A is a free spectral parameter. The linear systems (|34(l have two basis solutions ^ r± = {ipf, fpf), from which 
we build the so-called 'squared basis functions' 

f = -^t^2+Tpiri), g = i>ti>i, h = -^ 2 . (35) 

They satisfy the following linear systems 

f x = —ug — uh, g x — —2uf — 2i\g, h x — —2uf + 2i\h, (36) 

and 

ft = -iCg + iBh, g t = 2iBf + 2Ag, h t = -UCf - 2Ah, (37) 
and have the following integral 

/ 2 - gh = P(X) (38) 

independent of x and t. The periodic solutions are distinguished by the condition that P{X) be a polynomial in A and 
we shall see that one-phase solution corresponds to the sixth degree polynomial in even powers of A, 

3 

P(X) = Y[(X 2 - A 2 ) = A 6 — Sl A 4 + .s 2 A 2 - s 3 . (39) 

i=l 

Then /, g, h, satisfying (|36|) - (|38[) should be also polynomials in A, 

f = X 3 -f 1 X, g = iu(\- fii)(\- fin), h = -iu(X + m)(X + fj, 2 ), (40) 
where fij are new dependent variables. Substitution of (|40|l into (|38|l gives the conservation laws (sj are constants) 

2/i+u 2 = Si, f?+u 2 (n 2 1 +(4) = s 2 , u 2 ^ 2 ^ s 3 , (41) 
and substitution into (|36[1 and l|37|l yields the following important formulae 

u x = 2iu(ni + /j, 2 ), (42) 

u t = 2(2/i+« 2 )« !e . (43) 
^From (|43H and the first equation i|41[) we see that u depends only on the phase 

u = u(9), 6 = x + 2s 1 t, (44) 
and from (|42|l and the other equations (|41|l we find 

u 2 e = u A - 2siu 2 T + 4- 4s 2 = Q(u), (45) 
where the zeroes Ui of the polynomial Q{u) are related with the zeroes Xi of the polynomial P(X) by the formulae 

ui = ±(A 1 +A 2 + A 3 ), u 2 = ±(Ai - A 2 - A 3 ), u 3 = ±(-X 1 + X 2 - A 3 ), u 4 = ±(-Ai - A 2 + A 3 ). (46) 
If we order the zeroes A^ according to 

Ai > A 2 > A 3 > A 4 , (47) 
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then for the upper choice of the sign in 1|45[1 and (|46|l we have 

ui > u 2 > 113 > 1*4 (48) 

and u oscillates within the interval 

u 3 <u< it 2 , (49) 
where Q{u) > 0. For the lower choice of the sign in 1|45|) and (|46J) we have 

Ml < U2 < Us < 1*4 (50) 

and it oscillates within the interval 

u 2 < u < 1*3. (51) 

We are interested in wave trains against positive constant background which corresponds to the lower choice of sign 
in l|45l) and (|46l) . In this case Eq. (|45|l yields the periodic solution 



ta\ ("3 - ""1)112 - («3 - u 2 )uisn 2 (y/ (u 4 - u 2 )(u 3 - m) 9/2, m) 

u(9) = = , (52) 

u 3 - ui - {U3 - u 2 )sn 2 (y / (114 - u 2 ){u 3 - ux) 8/2, m) 



where 

_ (lt 3 - 1t2)("4 - ui) - Xl 



x + 2sit = x + 2{X\ + X\ + A|)t. (53) 



(m 4 - u 2 )(u 3 - m) A^-A§' 
At A2 = A3, when m = 1, the solution (|52|l transforms into soliton solution of the mKdV equation 

Us(0) = X 1 -- n 2 } Xl ~ 2 }} /v^^ T, e = x + 2{Xl+2Xl)t. (54) 

Ai - A 2 + 2A 2 cosh (2^/Xl - Xl 6) 

In a modulated wave the parameters A^ become slow functions of x and t. It is convenient to introduce new variables 



ri = A^ r 2 =A2, r 3 = X\, (55) 
so that Whitham equations can be written in the form (see, e.g. 0) 

+ = 0, «i=(l-A*)v, i = l,2,3, (56) 



dt dx \ diL 

where V = — 2{r\ + r 2 + ^3) is the phase velocity of the nonlinear wave l|52|l and 

K(m) r 3 - r 2 

L = - m = (57) 

V r 3 -r% r 3 - ri 

is the wavelength. 

Now our task is to consider the solution of the mKdV equation after the wave-breaking point. As it follows from 
(|20|1 . before this point in dispersionless approximation the evolution of the pulse obeys the Hopf equation 

Ut — 6u 2 u x = (58) 

with the well-known solution 

x + 6u 2 t = /(it 2 ), (59) 

where /(it 2 ) is determined by the initial condition. At the wave-breaking point, which will be assumed to be t = 0, 
the profile u 2 (x) has an inflexion point with vertical tangent line, 
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u 2 (x) 









t = o/ 


t>o\ 





FIG. 2: Formation of multi- valued solution of the mKdV equation in dispersionless limit J^Bj. The initial data correspond to 
2 - 2n3 (see Eq. 



the cubic curve a; = — (it 2 — it 2 .) 3 



Hence, in its vicinity we can represent l|59|l as 



x + 6u t 



> 2 



,2\3 



(60) 



where w& = u(xb,tb). Note that the mKdV equation is not Galileo invariant and therefore we cannot eliminate the 
constant parameter u 2 , in contrast to the case of a KdV equation (see, e.g. 

For t > solution HfiOjl becomes multi- valued function of x. Formation of this multi- valued region is shown in 
Fig. 2. For t > we cannot neglect dispersion and have to consider full mKdV equation. Due to effect of dispersion 
the multi- valued region is replaced by the region of fast oscillations which can be represented as a modulated periodic 
solution of the mKdV equation (|21(1 . We rewrite this solution [see Eqs. Ij52(l . (1531) ] in terms of the slowly varying 
functions ri(x, t), i — 1,2,3, 



u(x, t) = 



where 



(y/fa + Vn)(V^ +Vn- 05) + Wri - V^)(Vn + V^+ 7^)sn 2 (2Vr 3 - fj 9, m) 

sfn + \Ai - ( - V^2) sn2 ( 2 V r 3 - n 9, m) 

= x + 2stt = x + 2(n +r 2 + r 3 )t. 



r3 - r 2 



r3 - n 



(61) 



(62) 



and functions fi{x, t) are governed by the Whitham equations (|56|l . We have to find such solution of these equations 
that the region of oscillations matches at its end points corresponding to m = and m = 1 to the dispersionless 
solution (|60|l which we rewrite in the form 



x + 6rt = — (i — rb) 3 , r = u 2 , Tb = it 2 . 
This means that the solution of Eqs. I|56|) written in implicit form 

x- Vi(r)t = Wi(r), i = l,2, 3, 

must satisfy the boundary conditions 



W l\r 2 =r 3 



:r 3 =- 6r i! w 3| r2=ri = ~6r 3 ; 
-(r!-r b ) 3 , w 3 \ r2=ri =-(n -r b f. 



(63) 
(64) 
(65) 
(66) 



Then, as we shall see from the results, mean values of u will match at these boundaries to the solution of dispersionless 
mKdV equation. 

To find the solution l)64|l subject to the boundary conditions H65f) . l|66|l . we shall follow the method developed earlier 
for KdV equation (see, e.g. pj)- We look for Wi in the form similar to 



>V; = [l-^di)W. / = 1.2.3. 



(67) 
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and find that W satisfies the Euler-Poisson equation 

9 " w -wh-) 



(diW - djW) = 0, i^j. 



(68) 



For our aim it is enough to know a particular solution of this linear equation W = const/ y/P(r), where P(r) is a 
polynomial with zeroes r*j and it can be identified with polynomial (|39|) with taking into account Eqs. (|55|l . The series 
expansion of this solution in inverse powers of r, 



W = 



_4 r 3/2 



5 Q 3 
*1 



y/{r-r 1 )(r-r 2 )(r-r 3 ) n=Q - 
can be considered as generating function of sequence of solutions, 

WW = -2s 1: WW=2s 2 -lsi, =3s lS2 -2s 3 - - ri 

where si, s 2 , s 3 are the coefficients of the polynomial (|39() expressed in terms of r%, r 2 , r 3 : 

si = n + r 2 + r 3 , s 2 = nr 2 + rir 3 + r 2 r 3 , s 3 = rxr 2 r 3 . 
It is easy to find that the resulting velocities 



» 



L 



= ^-M di]w(n) ' 



1,2,3, 



have the following limiting values 



,(!) 



U llr 2 =r 3 ~ 6r l> W 3 



(1) 



-6r 3 



(69) 

(70) 
(71) 

(72) 
(73) 



15^2 



= — T?r 



„( 2 ) 



_ 2 r 



3' 



(74) 



(3) 



35 „3 



1 ■ 



Thus, we see that if we take 

Wi{r) 



W 



(r) + ±r b wf ] (r) - ±w<f> (r), i = 1, 2, 3, 



(3) 



35 „3 
X r 3- 



(3) , 



(75) 



(76) 



then formulae l|64|l satisfy all necessary conditions and define dependence of ri,r 2 ,r 3 on x and t in implicit form. In 
Fig. 3 we have shown the dependence of n, r 2 , r 3 on x at t = 0.5 and rb = 5. It is clearly seen that r 2 and r\ coalesce 
at the right boundary x+, where m = 1, and r 2 and r 3 coalesce at the left boundary where to = 0. Dispersionlcss 
solution is depicted by dashed line and r\ matches this solution at x_ and r 3 matches it at x + . 

Let us find the laws of motion x±(t) of the boundaries of the region of oscillations. At the right boundary we have 
the condition 



dx 



dx 
dr 2 



= 0, 



which yields the expression 



12 „2 
S5 r l 



35 



3 r 2 
70 r 3 



jin 



\r 3 r b + \rl, 



(77) 



(78) 



and substitution of this expression into Eqs. I|64() with 7*1 = r 2 gives the coordinate x+ expressed in terms of the 
Ricmann invariants r\ and r 3 : 



± [-32rf + 2r 3 3 - 8nr 3 (r 3 - 7r b ) - 8r 2 (4r 3 - 7r b ) - lr\r h 



35r, 3 l 



(79) 
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r 1 <r 2 <r 3 




FIG. 3: Dependence of the Riemann invariants ri,ri,rz on x at some fixed value of time. The plots are calculated according 
to formulae I64|l with r\> = 5 and t = 0.5. A dashed line represents the dispersionless solution which matches the Riemann 
invariants at the boundaries x± of the region of oscillations. 



On the other hand, this value of x+ must coincide with coordinate obtained from the dispersionless solution (|63[) with 
t equal to Eq. JTSJl, 



x+ = -6r 3 t - (r 3 - nf = ^ [-72r 2 r 3 + 26r| - 24nr 3 (r 3 - 7r b ) - 63r 2 r b - 35rg] 
Comparison of these two expressions for x+ yields the relation between r\ , r 2 and r 3 at m = 1 : 

(4ri + 3r 3 ) m=1 = 7r 6 . 
Substitution of r x obtained from this equation into Eq. (|78|) gives 



Hence 



and again with the use of Eq. (|81|l we obtain 



r 3 | m= i =r fc + |Vl5i 



r iL=i = r 2 | m =i = r 6 - ±\/l5l 



(80) 
(81) 
(82) 
(83) 
(84) 



These formulae give values of the Riemann invariants at the right boundary as functions of time t. Their substitution 
into (|79l) or l|8U|l yields the motion law of the right boundary 

(85) 



X+(t) = -6r b t + | J 1 1 3 / 2 . 



In a similar way at the left boundary x- the conditions 



dx 
dr 2 



dx 

i=0 dr 3 



= 



m— 



yield 



30' 1 15' 11 J 1 5' 3 1 3' 11 3 

which substitution into Eqs. (|64fl and I|rj3|) gives, respectively 



(86) 



(87) 



1 [-2rf - 32r| + 8rir 3 (4r 3 - 5r h ) + 40r 2 r b - r 3 b + r 2 (-8r 3 + 15r 6 )] 



and 



X- = \ [&r\ + r 2 (8r 3 - 25r b ) - 8nr 3 (3r 3 - 5r b ) - 5r 3 ] 



(89) 
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Their comparison yields the relation 

(?-i+4r 3 ) m=0 = 5r b (90) 
which permits us to eliminate r 3 from Eq. (|87J) to obtain 

t^^in-r,) 2 . (91) 

Hence 

ri| TO=0 =r 6 -2V3t (92) 

and again with the use of Eq. (|81|l we obtain 

r 2 | m=0 = r z \ m=0 = r b + \Vit. (93) 
Substitution of these formulae into Ij88(l or (|89|l yields the motion law of the left boundary 

x_ (t) = -Qr b t - 12V3 1 3/2 . (94) 

The plots of x+(t) and x-(t) are depicted in Fig. 4. To the right from x+(t) and to the left from x-(t) the wave 
is described by the dispersionless solution (|63|) . Between x+(t) and X-(t) we have the region of fast oscillations 
represented by Eq. I|62|l with ri{x, t), i = 1, 2, 3, given implicitly by Eqs. 1761) . The dependence of u on x at some fixed 
moment of time is shown in Fig. 5. It describes dissipationlcss shock wave connecting two smooth regions where we 
can neglect dispersion effects. At the right boundary the periodic wave tends to a sequence of separate dark soliton 
solutions of the mKdV equation and at the left boundary the amplitude of oscillations tends to zero. 




FIG. 4: Dependence of coordinates x± of boundaries of the region of oscillations on time for the mKdV equation case. 



B. Dissipationless shock wave in the KdV(2) equation I1M2H 



The theory of dissipationless shock wave for the KdV(2) equation is similar to that for the KdV and mKdV cases. 
Therefore we shall present here only its main points. 

The periodic solution of KdV(2) equation has the same form as in KdV equation case (see, e.g. |f|), 

u(x, t) = T2 + r3 — n — 2(r2 — 7'i)sn 2 (V r 3 — r i ®i m ) ; m — — -, (95) 

r3 — ri 

with phase velocity 

9 = x-Vt, V = 2s 2 -^sl, (96) 
corresponding to the second equation of the KdV hierarchy. 
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FIG. 5: The dissipationless shock wave for the mKdV equation. The parameters are equal to ri, = 5 and t = 0.5. The dashed 
line represents the squared roots of the Riemann invariants which match the smooth solution u(x, t) of the dispersionless 
equation at the boundaries of the region of oscillations. 

Now the periodic solution l|95|) is parameterized by the Riemann invariants r^, i = 1,2,3, rather than by their 
squared roots, as it was in mKdV equation case. Hence, the solution of the dispersionless equation 

u t = fu 2 u x (97) 

near the wave breaking point should be taken in the form 

x+^-u 2 = -{u-u b f. (98) 

In fact, this form is equivalent near the wave breaking point to the solution (|60|l . since {u 2 — u 2 ) ~ 2u b {u — u b ) and the 
constant factor can be scaled out. Formation of multi-valued region is illustrated in Fig. 6. After taking into account 
the dispersion effects it should be replaced by the dissipationless shock wave. 



U(X) 




X 



FIG. 6: Formation of multi- valued solution of the KdV(2) equation in dispersionless limit l|97|l . The initial data correspond to 
the cubic curve x = — (u — ut,) 3 (see Eq. 1981 ). 

Within the shock wave we have modulated periodic solution (|95|l where are slow functions of x and t and their 
evolution is governed by the Whitham equations Ij56|l with V defined by Eq. I|96|l . Their solution subject to the 
necessary boundary conditions can be found by the same method as was used in the preceding subsection. As a result 
we obtain 

x - wfh = ul + - lu bW f ) + ±w?\ (99) 

where are defined by formulae l|7(J|l - H72|l . Equations l|99|) define implicitly the dependence of the Riemann 

invariants r^, i = 1,2,3, on x and t. The resulting plots are shown in Fig. 7. At the right boundary x + we have 
soliton limit (to = 1) and at the left boundary a;_ we have a wave with vanishing modulation. 
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r 1 <r 2 < r 3 




FIG. 7: Dependence of the Riemann invariants r\,ri,ra on x at some fixed value of time for the KdV(2) equation case. The 
plots are calculated according to formulae l|99|l with Ub = —0.5 and t — 0.25. A dashed line represents the dispersionless 
solution which matches the Riemann invariants at the boundaries x± of the region of oscillations. 



t 




FIG. 8: Dependence of coordinates x± of boundaries of the region of oscillations on time for the KdV(2) equation case. 

The motion laws x±(t) can be found as in the preceding subsection, but the final formulae become now quite 
complicated and we shall not write them down. The corresponding plots are presented in Fig. 8. Again the region 
between X-(t) and x+(t) corresponds to expanding with time t dissipationless shock wave. It is illustrated in Fig. 9 
where the dependence u(x) at fixed moment t is shown. Now at one boundary bright solitons are formed and at the 
other boundary amplitude of oscillations tends to zero. 

IV. CONCLUSION 

We have shown that GDNLS equation with finite density boundary conditions can be reduced, depending on 
values of the parameters 77 and p, to several important continuous models — KdV, mKdV, KdV(2) and fifth-order 
KdV equations which describe different regimes of wave propagation in nonlinear lattice (Salerno model). The 
KdV(2) equation appears at such values of the parameters for which nonlinear and dispersive terms in in the KdV 
equation vanish, so that the main contribution in small amplitude long wave approximation is given by the third order 
nonlinear and fifth-order dispersion effects. This point correspond to the AL equation case, and since the AL equation 
is completely integrable and multi-scale method preserves complete integrability, we arrive at the second equation of 
the KdV hierarchy. This observation explains physical meaning of higher equations of integrable hierarchies — they 
give main contribution into wave dynamics if lower order effects vanish in small amplitude approximation of initial 
integrable equation. 

The evolution equations obtained as approximations to the GDNLS equation lead in dispersionless limit for general 
enough initial pulses to wave breaking so that taking into account small dispersion effects leads to formation of 
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t = 0.25 



u(x) 



-10 




-2 



-1 



2 



3 



X 



FIG. 9: The dissipationless shock wave for the KdV(2) equation. The parameters are equal to ut = —0.5 and t = 0.25. The 
dashed line represents the Riemann invariants which match the smooth solution u(x,t) of the dispersionless equation at the 
boundaries of the region of oscillations. 

dissipationless shock waves. We have developed the theory of these waves in framework of Whitham averaging 
method. Analytical expressions are obtained which describe their main characteristics — trailing and leading end 
points, amplitudes and wavelengths. 

The phenomena described in the present paper are not restricted by the GDNL, but are characteristic features of 
a large class of nonlinear Schrodinger lattices, which depend on one or more free parameters. 
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